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Abstract Solar observations in the infrared domain can bring important clues 
on the response of the low solar atmosphere to primary energy released during 
flares. At present the infrared continuum has been detected at 30 THz (10 ^m) 
in only a few flares. SOL2012-03-13 , which is one of these flares, has been 
presented and discussed in Kaufmann et al. (2013). No firm conclusions were 
drawn on the origin of the mid-infrared radiation. In this work we present a 
detailed multi-frequency analysis of the SOL2012-03-13 event, including observa¬ 
tions at radio millimeter and sub-millimeter wavelengths, in hard X-rays (HXR), 
gamma-rays (GR), Ha, and white-light. HXR/GR spectral analysis shows that 
SOL2012-03-13 is a GR line flare and allows estimating the numbers of and 
energy contents in electrons, protons and a particles produced during the flare. 
The energy spectrum of the electrons producing the HXR/GR continuum is 
consistent with a broken power-law with an energy break at ~ 800 keV. It 
is shown that the high-energy part (above ~ 800 keV) of this distribution is 
responsible for the high-frequency radio emission (> 20 GHz) detected during 
the flare. By comparing the 30 THz emission expected from semi-empirical and 
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Origin of the 30 THz emission during a solar flare 


time-independent models of the quiet and flare atmospheres, we find that most 
(^80%) of the observed 30 THz radiation can be attributed to thermal free- 
free emission of an op tically-thin source. Using the F2 flare atmospheric model 
( Machado et ai . 1980) this thin source is found to be at temperatures T ^ 8000 
K and is located well above the minimum temperature region. We argue that the 
chromospheric heating, which results in 80 % of the 30 THz excess radiation, can 
be due to energy deposition by non-thermal flare accelerated electrons, protons 
and a particles. The remaining 20% of the 30 THz excess emission is found to 
be radiated from an optically-thick atmospheric layer at T 5000 K, below the 
temperature minimum region, where direct heating by non-thermal particles is 
insufficient to account for the observed infrared radiation. 


Keywords: Radio Bursts, Microwave; X-Ray Bursts, Association with Flares; 
X-Ray Burst, Spectrum; Chromosphere, models; Heating, Chromospheric; Heat¬ 
ing, in Flares 


1. Introduction 


Solar flare continuum observations from the mid-infrared domain (a few tens of 
THz) to the millimeter-sub-millimeter radio domain provide in principle unique 
diagnostics of energy transport processes from the flare energy release region to 
the chromosphere_and_ofthejnost_energeticflm;e^ecelerated£articles_[e^Trot- 
tet and Klein, 20131 . and references therein'). Since year 2000. radio observations 
at 212 and 405 GHz have been routinely obtained by the Solar Submillimeter 
Telescope (SST) and more recently high-cadence imagi ng observations of both 


the quiet Sun and flares have been obtained at 30 THz (iKaufmann et ali 12008 
and references therein). The 2012 March 13 flare at ^ 17:22 UT ( S OL2012-03- 


13) was the hrst flare observed at 30 THz (|Kaufmann et a/.! . 120131 1. Since then 


othei_flnres^etected_a4.30_THzJiave^eenn;e2Oi;tedJnAheJiterature_(Kaufmann 
et aL, l2015l. and references therein). During SOL2012-03-13 the time evolution 
of the 30 THz impulsive emission was found to be similar to that of the radio 
(1-212 GHz), hard X-ray and white-light emission. The 30 THz emitting source 
was also co-spatial with a region of flare-enhanced white light continuum. No 
firm conclusion on the origin of the 30 THz inf r ared emission observed during 
SOL2012-03-13 was drawn bv IKaufmann et al\ ( 2013h . Nevertheless, they sug¬ 
gested that it is consistent with heating of the flaring atmosphere below the 
temperature minimum region which cannot be due to direct heating by hard 
X-ray emitting electrons, at leas t within the classic thick -target model. On the 
other hand, theoretical work bv lOhki and HndsonI (1975) raised the possibility 
that the infrared continuum at wavelengths shorter than 20 ^m may arise 
from an optically-thin source at temperatures in the range of 10"^-10® K. This 
latter possibility is in line with the results of time-dep endent simulations of th e 
chromospheric heating by electron beams performed bv lKasparova et 'ai,. (2009a). 

The main goal of this paper is to investigate what is the origin of the in¬ 
frared continuum detected at 30 THz during SOL2012-03-13. For that we 
have performed a detailed analysis of multi-frequency observations including 
measurements in radio millimeter and sub millimeter wavelengths, hard X-ray 
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(HXR) and gamma-ray (GR), Ha, and white-light (WL). The instrumentation 
used to collect these data is briefly described in Section [2] and an overview of 
the observations is presented in Section |31 Section 2] is devoted to the analysis of 
HXR/GR spectra from which numbers and energy contents of electrons, protons 
and a particles accelerated during SOL2012-03-13 are estimated. In Section [5] 
we discuss the relationship between HXR/GR and radio emitting electrons. The 
origin of the 30 THz emission is examined in Section [5] by using semi-empirical 
time-independent models of both the quiet and flare chromosphere. In Section [7] 
we provide an estimate of the energy deposited by flare accelerated electrons, 
protons and a particles in the atmospheric layers which radiate the 30 THz 
emission and we examine if it can account for the energy radiated at 30 THz. 
The final Section summarizes the main conclusions. 


2. Instrumentation 


Radio data at 45 GHz and 90 GHz were obtained by the radio polarimeters 
described in Valio et al. (2013). For the SOL2012-03-13 flare the uncertainty 
on flux densities is about 20%. Observations at 212 and 405 GHz have been 
provided by the SST. The atmospheric opacity was estimated to be 0.47 and 2.2 
Nepers, at 212 and 405 GHz respectively. At 212 GHz the flux density accuraw 
is estimated to be 20%. At 405 GHz the 1 a detection threshold was ~ 55 sfio 
and no significant emission was detected. 

In the HXR and GR domain, the SOL2012-03-13 flare has been detected by 
the Gamma-Ray Burst Monitor (GBM) on Fermi ( MeeEfan et o,l. . 2009l l. GBM 
is composed of twelve sodium-iodide (Na I) and two bismuth germanium oxide 
(BGO) detectors. In the present analysis we have used the most sunward Na 
I and BGO detectors. The Na I data consist in 128 energy channels covering 
the energy range 4 keV-2000 keV and the BGO data also consist of 128 energy 
channels but in the 113-50,000 keV range. 


The 30 THz instrumental s etup at El Leoncito is reported in iMarcon et al. 


1)20081 1 and iKaufmann et a,l\ ()2008l l . For the present solar flare the data cali¬ 
bration and the determination o f the excess flux density relative to quiet Sun 
have been described in details in IKaufmann et all (2013|). The accuracy of flux 


density measurements has been estimated to be ~ 25%. The diffraction limit 
of the instrument is ^ 15" and the source position uncertainty is cr 3" as 
obtained from fitting of the 30 THz limb. The bandwidth of the instrument is 
^17 THz and in the following, we assume that the 30 THz emission does not 
depend on frequency within this bandwidth. 

Ha observation s were obtained wit h the Fl-Alpha Solar Telescope for Ar¬ 
gentina (HASTA; iBaeala et al. . 19991) . HASTA provides solar full disk images 
with a time resolution of 5 seconds in flare mode and a spatial resolution of 2". 
Excess emission due to the flare has been normalized to the quiet Sun emission. 

The above data set has been complemented by radio observations at 1.415, 
2.695, 4.995 and 8.8 GHz from the Sagamore Hill-station of the US Air Force 


Usfu = . m-2 ■ Hz- 
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Figure 1. Temporal evolution of the flare SOL2012-03-13 in various wavelength domains 
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Figure 2. Projections of Ha (17:22:53 UT, red contours), WL (17:23:00 UT, yellow contours) 
and 30 THz (17:22:54 UT, green circle) sources on a SDO/HMI vector magnetogram obtained 
at 17:13:15 UT 


Radio Solar Telescope Network (RSTNj^ (iGuidicj . 1979f) . During SOL2012- 
03-13 there were no reliable 15.4 GHz measurements. We have also used WL 
images from the Helioseismic and Magnetic Imager (HML [Scherrer et al. . 2012 ) 
on board the Solar Dynamic Observatory (SDO) obtained every 45 s with a 
spatial resolution better than ^ 1 " H. 
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3. Observation overview 


The flare SOL2012-03-13 occurred in NOAA active region AR 11429 located 
at N18W64 on 2012 March 13 at 22:00 UT. It was associated with a M7.9 soft 
X-ray (SXR) event detected by the Geostationary Operational Environmental 
Satellite {GOES) and with a Ha flare of importance IB. From top to bottom 
Fig. m displays the time profile of the SXR (1-8 A) from GOE^ radio (45 and 
212 GHz), Ha, flare excess mid-infrared (30 THz), WL, and 216-318 keV HXR 
emissions detected during SOL201 2-03-13. Except Ha me asurements, these 
observations have been presented in iKanfmann et al\ (|2013r i . The HXR light 
curve consists in an impulsive burst which started at ^ 17:20:53 UT with a 
maximum at 17:22:54 UT and lasted for about ten minutes. This impulsive 
peak is seen from radio to optical wavelengths, i.e. from the corona down to 
the low chromosphere. At radio wavelengths, there are two peaks of apparently 
similar flux densities. There is a change of the decay slope of the HXR emission 
which corresponds to the second radio peak. The Ha, 30 THz and WL emissions 
show a broad peak which encompasses the two radio and HXR peaks. In addition 
to this impulsive phase, there is a long-lasting and slowly varying component seen 
at all wavelengths except in HXR and WL. For example at 45 GHz, although it 
is weak (~ 40-50 sfu), the gradual emission is still visible after 19:00 UT. 

Figure [2] shows the projections of the Ha, WL and 30 THz sources, obtained 
around the maximum of the impulsive peak, on a vector magnetogram recorded 
by SDO/HMl at 17:13:15 UT. The flare enhanced Ha emission arises from four 
sources (red contours) marked Si, S 2 , S 3 and S 4 . Si and S 2 are close to the two 
WL emitting regions, WLi and WL 2 , shown by yellow contours (see also Figures 
3 and 10 in iKaufmann et al. . 2013). The larger extents of WLi and WL 2 are 
respectively 10" and 2". The projections of WLi and WL 2 and Si and S 2 are 
located on opposite magnetic field polarities and are connected by a bright loop 
shown at 211 A in Figure 5 of iKaufmann et a,l\ (2013). At 30 THz only the source 
corresponding to SI and to WLi, the stronger WL emitting region, is detected 
within the sensitivity of the present setup of the 30 THz instrument. Since the 
30 THz source is not resolved by the instrument, it is marked by a green circle 
whose diameter is the diffraction limit of the instrument (15") in both directions. 
Figure[2]thus indicates that the 30 THz source and the WL barycenter emission, 
which is close to maximum of WLi, are located close to each other. However, 
the uncertainty on the 30 THz source position (3 cr ~ 9") prevents to draw 
any conclusion on the relative height of both sources, for relative heights lower 
than few thousands of km. While the time profiles of Si and S 2 show a clear 
impulsive phase counterpart (see Figure [T|), the remote sources S 3 and S 4 only 
show a slowly varying emission that will not be considered in the following. 


^http: //www. ngdc.noaa.gov/stp/space-weather/solar-data/solar-features/solar-radio/ 
rstn-1-second/ 

®http: //hmi.stanford.edu/Description/HM LOverview.pdf 

^provided by NASA/GSFC at http://umbra.nascom.nasa.gov/goes/fits/ 
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Energy [keV] 

Figure 3. HXR/GR spectrum around the maximum of SOL2012-03-13 obtained by combin¬ 
ing BGO and Na I measurements of the GBM instrument on Fermi (error bars). The green 
and red curves show the model spectrun fitted to GBM observations (see text). 


4. Electron and ion numbers and energy contents around the 
maximum of SOL2012-03-13 

Figure [3] shows the mean excess count rate spectrum in the 20-10,000 keV range 
recorded by the most sunward Nal (black error bars) and the most sunward 
BGO (blue error bars) detectors, for a 80 s time interval around the peak of 
the HXR/GR emission. In the 36-1000 keV range, the observed Nal count rate 
spectrum is well represented by the expected count rate spectrum shown by 
the green line. This latter spectrum has been obtained by convolving a triple 
power-law photon spectrum with the Nal detector calibration matrix. The triple 
power-law spectrum represents the electron bremsstrahlung spectrum in the ^ 
35-1000 keV range. It is defined by the normalization factor at 50 keV, the 
power-law indices 71 , 72 and 73 and the break energies Ebri and Ebr 2 displayed 
in Table [H In the 216-10,000 keV range, the BGO count rate spectrum is well 
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Table 1. Parameters of the photon spectrum fitted to the Nal measurements in 
the 35-1000 keV 


photons/s/cm^/keV 71 

at 50 keV 

Ebri 

keV 

72 

Ebr 2 

keV 

73 

4.7 ± 1.5 

Table 2. Parameters 
200-10,000 keV 

3.3ib 0.2 94 zb 10 3.8 ±0.2 282 ± 20 2.4 ± 0.2 

of the photon spectrum fitted to the BGO measurements in the 

photons/s/cm^ /keV 
at 50 keV 

Pi 

Ebra 

keV 

r2 

2.2 MeV line 
photons / s 

narrow lines 
photons/s 

5.7 ± 1.9 

3.6 ± 0.2 

388 ± 50 

2.3 ± 0.2 

0.054 ± 0.018 

0.11 ± 0.04 

electrons / s/keV 

5i 

Ebre 

^2 



at 50 keV 


[keV] 




5.6 X 10^3 

4.7 

800 

3.5 




represented by the expected count rate spectrum (red line) obtained for a photon 
spectrum which is the sum of a broken power-law, of the 2.2 MeV neutron capture 
line and of a template of narrow nucl ear lines. We used the nuclear line template, 
included as standard within OSPEX ([Schwartz et oi.l . l2002h . calculated for a flare 
at an heliocentric angle of 60° by assuming a downward isotropic distribution of 
ions with a power-law energy distribution of index 4 and an a/p ratio of 0.22. 
The broken power-law represents the electron bremsstrahlung spectrum in the 
200-10,000 keV range. It is defined by the normalization factor at 50 keV, the 
power-law indices Pi and P 2 and the break energy Ebra (see Table [2|) . For the 
fitting procedure we have used the OSPEX package of Solar Soft. 

Figure [3] shows that: (i) the electron bremsstrahlung continuum extends up 
to at least 10,000 keV and (ii) there is clear 2.2 MeV line emission and weak 
narrow GR line emission indicating that ^ 1-100 MeV protons and ^ 1-100 
MeV/nucleon ions have been accelerated during SOL2012-03-13. The photon 
fluxes derived from both Na I and BGO data agree well (within better than 
25%) between 200-700 keV. However, because the responses of Nal and BGO 
detectors are quite different, a given incident photon spectrum will produce 
different count spectra in both types of detectors as shown in Figure^] Above ^ 
400 keV, the l-tr uncertainty is smaller for BGO than for Na I. In the following 
we will only consider the photon spectrum derived from BGO data since we 
are interested in microwave emitting electrons and in GR line radiation. In 
order to derive the spectrum of electrons from this photon spectrum we further 
assume that this electron spectrum is also a broken power law. The electron 
spectrum is then derived by using an electron-to-thick-target bremsstrahlung 
code including both relativistic effects and electron-electron bremsstrahlung 


®N. Vilmer, private communication 


SOLA: trottet_et_al_20120313_R2.tex; 22 September 2015; 0:32; p. 9 













Trottet et al. 


and by varying its parameters until a good agreement with the photon spec¬ 
trum is reached. The obtained parameters are displayed in Table [21 The mean 
flux of > 50 keV electrons between 17:22:10 and 17:23:30 UT is found to be 
Fe(> 50 keV) ^ 8 X 10^'^ electrons s“^ corresponding to a mean energy flux 
$e(> 50 keV) ~ 9 X 10^^ ergs . The nuclear line template is normalized 
so that an amplitude of 1 corresponds to 8.5946 x 10^® protons with energies 
above 30 MeV. Here the amplitude of the nuclear line template is 0.11 ± 0.04 
s“^ (see Table|2|). It is well documented that ions with energies down to almost 1 
MeV/nucleon contribute to the GR line emission. For the assumed energy power- 
law index of 4 we then get mean proton and a fluxes at the top of the atmosphere: 
Fp(> 1 MeV) ^ 2.6 x 10^^ protons s“^, Fa(> 1 MeV/nucl) 5.6 x 10^^ a s“^. 
These correspond to energy fluxes of respectively 4>p(> 1 MeV) = 6.2 x 10^^ 
ergs s“^ and $„(> 1 MeV) ^ 5.4 x 10^^ ergs s“^. Since the uncertainty on the 
determination of photon fluxes is about 30%, the above values of particle and 
energy fluxes suffer from the same uncertainty. The broad GR line emission has 
not been taken into account to get the above estimates because, due to the poor 
statistics, the observed GR line spectrum is not well defined. Nevertheless, we 
have checked that the inclusion of broad lines lead to estimates of proton and a 
particle fluxes and energy fluxes that stay within the 30% uncertainty. It should 
be noted that $p ^ and that <l’e(> 50 keV) ~ $p(> 1 MeV) + $a(> 1 MeV). 


5. Relationship between radio and HXR/GR emitting electrons 
during the impulsive phase 

Figure HI displays the averaged radio spectrum (vertical error bars) observed 
during the same time interval (17:22:10-17:23:30 UT) as the HXR/GR spectrum 
analyzed in Section [d] Its shape is reminiscent of a gyrosynchrotron spectrum 
with a turnover frequency somewhere between 10 and 30 GHz. It is generally 
admitted that the 1-200 GHz (cm-mm) radio emission and the HXR/GR contin¬ 
uum are both r adiated by a closely related population of electrons accel erated in 
the corona le.o. iBastian. Benz, and Gary . 19981 : Pick and Vilmeil . l2008f l. Briefly 
stated: (i) the gyrosynchrotron emission is radiated by the instantaneous pop¬ 
ulation of electrons present in the coronal portion of magnetic loops connected 
to the acceleration region, while (ii) the HXR/GR bremsstrahlung continuum is 
produced by thick-target interactions of precipitating electrons at the chromo¬ 
spheric foot points of these loops. It is also well documented that for magnetic 
field of a few hundred Gauss the cm-mm emission is radiated b y electrons with 
energies ranging from a few hundreds of keV t o a few MeV {e.g. iRamatv . 19691 : 


Pick. Klein, and Trottet . 1990l : Ramatv et al . 1994). As a first approximation 


we thus assu me that the instanta neous spectrum of radio emitting electrons is 
given by {e.g. Trottet et all 1998h : Nr(E) « Fx(E) Ttrap electrons keV“^ where 
Ttrap is the time spent by the accelerated electrons in the cm-mm radio emitting 
region and Fx(E) is the electron flux spectrum displayed in Table [2] Without 
multi-frequency cm-mm imaging observations it is not possible to constrain radio 
source models with inhomogeneous magnetic field B and ambient density Namb- 
In the following we thus compute the gyrosynchrotron emission produced by 
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Figure 4. Radio spectrum: in both panels the vertical error bars show the radio spectrum 
observed in the 1-212 GHz frequency range around the maximum of the impulsive phase of 
SOL2012-03-13: At 405 GHz the arrow indicates the upper limit of the flux density. The 
gyrosynchrotron spectra computed for B = 700 G and Ttrap = 2 s , Emax = 10, 000 keV and 
Emin = 30 keV (solid line), 400 keV (dotted line) and 800 keV (dashed line) are over plotted 
in the upper panel while results for Emin = 30 keV and Emax = 10,000 keV (solid line), 7,000 
keV (dotted line) and 4,000 keV (dashed line) are over plotted in the lower panel (see text). 
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Nfi(E) in a radio emitting region of constant Namb and uniform B. Such a 
simple model is not able to account for the < 10 GHz optically-thick radio 
emission, but it allows one, in principle, to reasonably match the optically-thin 
part of the radio spectrum observed at 45, 90 and 212 GHz during the present 
flare. The gyrosynchrotron emission has then been computed for B = 700 G, 


L trap 


= 2 s and Namb = 5 x 


10^° cm The diameter and thickness of the 


radio source normal to and along the line of sight have been set to 14,000 and 
2000 km respectively and the angle bet ween B and th e line of sight to 45° . 
The gyrosynchrotron code used is that by Ramatvl ( 1969l l: Ramatv et al. (1994) 
adapted to take into account a broken power-law electron spectrum. 

Figure 0] (upper panel) shows the results for an electron population with a 
high energy cutoff Emax = 10,000 keV and different low-energy cutoffs Emin = 
30 keV (solid line), 400 keV (dotted line) and 800 keV (dashed line). In Fig. 0] 
(lower panel), Emin = 30 keV and Emax = 10,000 keV (solid line), 7,000 keV 
(dotted line) and 4,000 keV (dashed line). Figure0]shows that the optically-thin 
part of the observed radio spectrum is emitted by electrons in the ^ 800 to 
7,000-10,000 keV energy range, that is above the break energy of the electron 
spectrum deduc ed from the HXR/GR photon spectrum. This is consistent with 
earlier findings ( Trottet et al. . 1998 . 2000l . 2008 ). B and Ttrap, which fix the value 
of Nr, are not independent parameters in the data fitting process. Indeed, the 
same flux densities can be obtained for larger B and lower Ttrap and vice versa. 
B ~ 700 G and Ttrap 2 s appears to be a reasonable set of values. On one 
hand, a substantial increase of B would lead to a too high turnover frequency. 
On the other hand, a decrease of B would imply an increase of Nr. Since there 
is no physical reason to have Emin larger than a few tens of keV, a too low value 
of B would lead to an unreasonably high density of non-thermal electrons in the 
radio source compared to Namb- 


6. Origin of the 30 THz emission 


In the following we assume that both the quiet and flare 30 THz emission is 
thermal emission from the chromosphere. The aim of this section is to search from 
which atmospheric layer(s) th e 30 THz emission obse rved during the SOL2012- 
03-13 flare (see Section 0] and iKanfmann et al. . 2013 ) arises. As we do not know 
the actual flaring atmosphere, to achieve this, a first approach employs time- 
independent models of both the flaring and quiet lower atmospheres of the Sun. 
Such semi-empirical models are based on hydrostatic and statistical equilibria. 
For a given temperature structure the full non-LTE radiative transfer is solved 
numerically in order to obtain the synthetic flare spectrum. This synthetic spec¬ 
trum is compared to observed spectral features (both lines and continuum) and 
the temperature structure is adjusted in order to get a good agreement with 
observations. 

In order to compare expectations from flare models with quiet Sun conditions 
we have used the C7 model for the quiet solar atmosphere (|Avrett and Looser . 


20081 ). The flare atmospheric models considere d here after are: the FLA model 
constructed by iMauas. Machado, and AvrettJ ( 1990l ) to account for the WL 
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Table 3. Height ( hmax), opacity (tsq) at the maximum of CF30. 
The altitudinal range Ah over which CF is larger than half of its 
maximum value. Tt,3o is the derived brightness temperature at 30 
THz. ATb3o is the increase of T[,3o with respect to C7 and AS30 
is the excess flux density with respect to C7 for a 30 THz emitting 
source of 10” diameter. 


Model 

hmax 

[km] 

T"30 

Ah 

[km] 

Tb30 

[K] 

2^Ti,3o 

[K] 

AS30 

[sfu] 

C7 

180 

1.1 

130-270 

4700 



FLA 

180 

1.1 

130-270 

4900 

200 

1000 

FI 

180 

1.4 

140-300 

5000 

300 

1500 

F2 

1070 

0.2 

960-1100 

7100 

2400 

12,300 


continuum emission observed during the SOL1983-06-15 flare, and the FI and 
F2 modelsre 2 resentatrve_ofres 2 ectivel^_weahandmedium;Size^are^Machado 


et al, 1980l l. We have restricted our analysis to these models because they are 
well documented in the literature and because SOL2012-03-13 is a medium size 
flare. We take their te mperature height profiles as inputs and use the PAKAL 
radiative transfer code ( De la Luz et al . , 20101 : IDc la Lu z. Lara, and Raulinl. 2011 ), 
which is similar to the PANDORA code used by Avrett and Loeserl ( 20031 '). in 
order to: (i) compute the densities of electrons and of the most important ion 
species, including H~, as a function of height h and, (ii) determine the absorption 
coefficient K 3 o(h) and the opacity ranfh) for a 30 TH z emitting region at a 
longitude of 60°. Following iHeinzel and AvrettI ( 2012h . we then compute the 
contribution function CF 3 o(h) as a function of h, given by: 


CF3o(h) = r;3o(h)e-"™(‘') , ry3o(h) = K3o(h)B3o(T(h)) 


( 1 ) 


where T(h) is the local temperature and B 3 o(T(h)) the Planck source function. 
CF 3 o(h) allows us to identify, for each model, the atmospheric layers which 
contribute most to the 30 THz emission. For a given model the brightness 
temperature at 30 THz is given by: 


Tb30 


Te "^^drso. 


( 2 ) 


For each atmospheric model Table 3 displays hmax and r 3 o at the maximum 
of CF 30 , the altitude range Ah at half of the CF 30 maximum, the expected 
brightness temperature Tbso, the excess brightness temperature ATbso and the 
corresponding excess flux density AS 30 (for a 10 " source) with respect to values 
obtained for C7. The dashed lines in Figure [S] show T 30 (top) and T (bottom) as 
a function of h for C7 (black) and F2 (red). The corresponding height prohles 
of CF 3 o(h) (solid lines) are over plotted in both panels. The 30 THz emission 
from C7, FLA and FI arises from an atmospheric layer LC7i centered at hmax ^ 
180 km, where r 3 o ~ 1, which extends from ~ 130-140 to 270-300 km . Thus, for 
these models, the 30 THz radiation arises from an optically thick layer located 
below the temperature minimum region at temperatures around 5000 K. For the 
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Figure 5. Top: Opacity rso (dashed lines) and contribution function CF (solid lines) as a 
function of height above the photosphere. Red curves represent F2 model results, while black 
curves are for C7. Bottom: Atmospheric local temperature (dashed lines) and CF (solid lines) 
as a function of height. 


quiet 


et aL, Il99li ). FLA and FI cannot explain the 30 THz emission observed during 
SOL2012-03-13 because they lead to AS 3 o~ 1000-1500 sfu, which is about one 
order of magnitude lower than the observed value 12,000 sfu at the maximum 
of the 30 THz burst. On the contrary, F2 leads to AS 30 12,300 sfu, which is 
in agreement with the observed value. In F2, the maximum of CF 30 occurs at 
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hmax 1070 km where rao < 1 so that the 30 THz radiation arises mostly (~ 
80%) from an atmospheric layer LF2, Ahh ^ 960-1100 km, located above the 
temperature minimum region at temperatures around 8000 K. Figure [S] shows 
that for F2, in addition to the LF2 layer, there is an optically thick 30 THz 
emitting layer, slightly below the temperature minimum, in the altitude range 
Ahi ~ 130-270 km. This latter source contributes moderately (~ 20%) to the 
total emission. Its temperature is a few tens of kelvin larger than that of the 
quiet Sun so that it represents only less than 1% of the quiet Sun emissio n. This 
is in line with earlier theoretical expectations bv lOhki and Hudson (1975) which 
indicate that the contrast between the quiet and flare infrared emission in the 
10-40 THz range is about ten times higher for an optically thin source at about 
10"^ K than for an optically thick source around 5000 K. This is also consistent 
with results of time-dependent numerical simulations ('see Figure 3 in Kasparova 
et aZ.. l2009ah . The F2 flare model appears thus as a plausible atmospheric model 
to account for the 30 THz emission radiated during SOL2012-03-13. 


7. 30 THz chromospheric response to flare accelerated particles 


The 30 THz time profile (see Figure [1} exhibits an impulsive emission, which 
encompasses the two radio and HXR peaks, superimposed on a long-lasting 
and slowly varying component (see Section [3|) . This suggests that the 30 THz 
impulsive burst results, directly or indirectly, from chromospheric heating due to 
energy deposition by non-thermal particles (electrons, protons and a particles) 
while the slowly varying emission is generated by some slower process like, 
conduction fronts. In the following we focus on the chromospheric heating by 
particles accelerated during SOL2012-03-13. A proper treatment of this problem, 
which is beyond the scope of the present study, should account for the temporal 
evolu tion of the flaring atmosp here due to time dependent heating by particles 


as m 


Kasparova et al\ (|2009al lbfl. Here, we simply consider C7 as the initial 
atmosphere model at to 17:20:45 UT (beginning of the HXR and 30 THz 
radiation) and F2 as the atmosphere model that is reached at ti 17:24:50 UT 
(maximum of the observed 30 THz radiation). For a given atmosphere model we 
have calculated the energy deposition per centimeter by non-thermal particles 
Q{N) (erg cm“^) as a function of the hydrogen column depth N (see Eq.[5]in the 
Appendix). In the following we use height h to parametrize altitude above the 
photosphere in a given atmosphere, considering 0(/i) rather than Q{N). These 
computatio ns are present ed in detail in the Appendix. They extend previous 
works {e.g. Emsliel . 1978h by considering that the degree of ionization varies 


with N and retaining the energy-dependence of the Coulomb logarithm, which 
is particularly important for ions. 

As a first approach, we have computed the total energy deposited by particles 
accelerated between to and ti as a function of h in C7. For that we have inte¬ 
grated the BGO count rate spectrum spectrum from to to ti and estimated the 
mean flux of electrons, protons and a particles in the same way as in Section 0] 
and multiplied by ti — to = 245 s . This leads to J'-oe(> 50 keV) ^ 10^^, 
J'op(> 30 MeV) ~ 8.5 x lO^o and J'oa(> 30 MeV/nucl) ~ 1.9 x lO^o for re¬ 
spectively the total numbers of > 50 keV electrons, > 30 MeV protons and 
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> 30 MeV/nucl a particles accelerated between to and ti. We then use these 
normalizations in Eq. (|S]) to calculate the energy deposited per cm due to elec¬ 
trons, protons and ions, and sum over these to obtain the total energy deposited 
by particles per cm E*°*(h). Figure [S] shows E*°*(h) as a function of h in the 
pre-flare atmosphere C7 for the ensemble of accelerated particles as well as the 
contributions from electrons, protons and a particles separately. As stated in 
Section [5] most of the 30 THz emission arises from the LF2 layer marked by red 
vertical dashed lines in Figured The altitude range of the C7 layer LC7 (750-820 
km), which covers the same range of mass column densities (3.5 x 10“^-6.2 x 10“^ 
g cm“^) as LF2, is marked by the black dashed lines in the Figure |6l Such a 
range of colum^densit^js_consisten^^it^^ha£^ndicate^^^Machado, Emslie, 
and Avrett Tl989^ . As a rough approximation we assume that the total energy 
deposited by particles between to and ti in LC7, Ep£’^, may account for the 30 
THz emission from LF2. By integrating E*°*(h) over the LC7 altitude range we 
get Ep^/ ~ 3 X 10^® erg. The total amount of energy radiated at 30 THz between 
to and ti is given by: 


7 

•/to 


E“‘ = 10-1^ 47r R^Aiy / AS 3 o(t)dt - 10^*^ erg. 


( 3 ) 


where R=1 AU=1.5 x 10^^ m, Aiy is the observed frequency bandwidth in Hz 
(see Section m and AS 3 o(t) the excess flux density at 30 THz in sfu. Ep^^ is 
about thirty times larger than Egg*. Thus, direct heating by non-thermal particles 
appears as a possible mechanism to account for the ^ 80% of the impulsive 30 
THz emission which arises from LF2. This latter statement remains valid even if 
we included a range of initial pitch angles, mirroring and pitch angle scattering 
in the energy deposition calculation (see Appendix). On the other hand, the 
energy 6x10^® erg deposited by particles in LC7i, which corresponds to about 
the same mass column density and altitude ranges as the equivalent layer in 
F2 (see Section IH]) is too small to account for the remaining 20% of the 30 
THz radiation. The heating in the LC7i layer may thus be due to some indirect 
mechanism, suc h as radiative backwarming of the deep chromosphere from the 
LF2 layer (e.ff., Metcalf. Canfield, and Saba . 19901 : Kerr and Fletcher . 2014 1 or 
radiative coupling of the upper chromo sphere and temperature minimim regions 
( Machado. Emslie. and Avrettl 1989h . It should be noted that a harder ion 
distribution {e.g. S = 2) would deposit a greater fraction of its total energy 
content in the deeper atmosphere. To this extent the calculations here, with a 
fixed value of 5 = 4 should be regarded as illustrative. The statistical quality of 
the 7 -ray spectrum for this event does not justify more detailed investigation, 
however. 

Figure [T] shows that the evolution of the Ha emission from the Si and S 2 
kernels (see Figure^]) is similar to that of the 30 THz excess flux density. Because 
the Ha emission arises from a ^ 10^ K plasma, the altitudes of Si and S 2 should 
be close to that of the 30 THz source at ^ 8000 K. This suggests that direct 
heating by particles is also responsible for the impulsive Ha radiation while, like 
for the 30 THz emission, the slowly-varying, long-lasting component is due to 
some other process. A detailed and quantitative analysis of the WL emission 
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Figure 6. Energy deposition per unit of height in C7 (E^°^) by flare accelerated particles as a 
function of height. The dashed vertical lines in red show the F2 atmospheric layer (LF2) from 
which ~ 80% of the 30 THz emission arises. The black vertical lines show the C7 layer (LC7) 
which covers the same mass column density as the LF2 layer. 


is beyond the scope of this study. Nevertheless, as remarked in Section [3l the 
WL time profile exhibits only the impulsive emission and not the slowly-varying, 
long-lasting component. This suggests that the mechanism responsible for the 
30 THz and Ha slowly-varying, long-lasting component is not efficient enough 
to produce WL emission whether it arises from direct particle heating from a 
LF2-like layer or from indirect heating of the deep chromosphere (LC7i layer). 
The former hypothesis is consistent with th e heights of WL sour ces, ^ 800-1000 
km above the photosphere estimated by e.g. Krucker et al. 1 2015l l for some flares. 


8. Summary 


The SOL2012-03-13 flare is the first event that have been detected in the mid- 


infrared domain at 30 THz (IKaufmann et a/.l . 1201311 . The combined analysis of 


the 30 THz data with HXR/GR, WL, Ha and 1-400 GHz radio observations 
has allowed us to estimate the numbers of electrons, protons and a particles 
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accelerated during the flare and to discuss the origin of the 30 THz emission. 

The main findings can be summarized as follows: 

The HXR/GR observations of the Gamma-ray Burst Monitor on FERMI 
have revealed that SOL2012-03-13 is a gamma-ray line flare detected up to 
at least 10 MeV. The numbers and energy contents of HXR/GR emitting 
electrons, protons and a particles have been estimated. The energy content 
in protons and a particles with energies above 1 MeV/nucl. is found to be 
similar to the energy content in > 50 keV electrons. This indicates that, in 
addition to electrons, protons and a particles may substantially contribute 
to the heating of the chromosphere during flares such as SOL2012-03-13. 
The observed 1-200 GHz radio spectrum is reminiscent of gyrosynchrotron 
radiation. The optically-thin part of the radio spectrum (> 10-20 GHz) is 
radiated by HXR/GR emitting electrons with energies greater than ^ 800 
keV. The spectrum of these high energy electrons is harder than that of 
lower energy electrons. 

Assuming that the 30 THz emission is thermal emission from the chromo¬ 
sphere we have used time independent semi-empirical models of both the 
quiet and flaring atmospheres in order to identify from which atmospheric 
layers the 30 THz emission arises. In agreement with earlier works we find 
that the quiet Sun 30 THz radiation is emitted by an optically thick layer, 
at temperatures around 5000 K, below the temperature minimum region. 
During the flare, the measured flux densit y excess at 30 THz is c onsistent 
with that expected from the F2 model bv iMachado et (il. 1 1980h which is 
representative of medium size flares. In this model, about 80% of the 30 THz 
emission arises from an atmospheric layer LF2 (~ 960-1100 km above the 
photosphere) at temperatures around 8000 K well above the temperature 
minimum region. At 30 THz, this layer is not optically thick (r < 1). The 
remaining 20% of the 30 THz emission arises from an optically-thick layer, 
slightly below the temperature minimum region in F2, which covers similar 
altitude and mass column density ranges to those of the emitting layer in 
the quiet Sun. 

The energy deposited by electrons, protons and a particles in the quiet 
Sun layer which covers the same range of mass column density (3.5-6.2 
xI0“^ g cm“^) has been calculated by taking into account that the degree 
of ionization varies with height and by retaining the energy-dependence 
of the Goulomb logarithm. This energy is found to be about thirty times 
greater than that radiated at 30 THz. Thus energy deposition by accelerated 
particles is a plausible mechanism responsible for the heating of the LF2 
layer which gives rise to 80% of the observed 30 THz emission during 
the impulsive phase of SOL2012-03-13. On the other hand, the energy 
deposited by particles in the deep chromosphere is too small to account 
for the remaining 20% of the 30 THz flux density radiated from below the 
temperature minimum region in F2. A possible, but not unique, scenario is 
that the impulsive WL emission and 20% of the 30 THz radiation arises from 
an optically-thick layer, below the temperature minimum region, heated by 
e.ff. backwarming from the upper chromosphere where most of the energy 
transported by accelerated particles is deposited. 
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The 30 THz and Ha exhibit a long-lasting and slowly varying component 
which is not observed in the WL and HXR/GR domains. This indicates 
that the chromospheric heating giving rise to this component is due to 
a different mechanism than energy deposition by accelerated particles and 
that the efficiency of this mechanism is not sufficient to produce long-lasting 
WL emission. 


In summary, using time-independent semi-empirical models of the quiet and 
flare atmosphere, we have shown that thermal emission from the chromosphere 
heated by accelerated particles provides a plausible quantitative interpretation of 
the 30 THz emission in gross agreement with that observed during the impulsive 
phase of SOL2012-03-13. Time dependent n umerical models of energy deposition 
by accelerated particle such as developed by Kasparova et all ( 2009al l are needed 
to confirm our conclusion. Moreover, if the 30 THz emission arises from a quasi 
optically-thin chromospheric layer like LF2, the emission spectrum between ^ 
10-100 THz is expected to be rath er flat or even slightly de creasing with in¬ 
creasing frequencies (see Figure 2 in Ohki and Hudson . 197fjl . Observations at 
frequencies above 30 THz such as those recently obtained from the McMath- 


Pierce Solar Facility at Kitt Peak (jPenn el o,ll l2015f) should allow us to test 


the validity of such an expectation, and thus of the thermal interpretation of 
the infrared emission from flares proposed in this work. In the near future, the 
combination of McMath observations and measurements at 3 and 7 GHz by 


Solar- T (|Kaufmann et ali . 120141 1 should allow to complete our understanding of 


the chromospheric response to flare energy release. 
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Appendix 


In order to compute the energy deposited by accelerated electrons, protons and a 
particles in the chromosphere as a function of depth, we inject particles with an 
energy spectrum Fq{E) (MeV“^) at the top of the atmosphere. We parametrize 
position in the atmosphere using the hydroge n column dep th N. The energy 
deposited (erg cm“^) at depth N is (following lEmslie . 1978 ) 


0(iV) 



Fo{Eo) 

v{E) 






(4) 


SOLA: trottet_et_al_20120313_R2.tex; 22 September 2015; 0:32; p. 19 
































Trottet et al. 


Here Emin{N) is the minimum energy of particle that can reach depth N. In 
Equation (|1]) E should be understood as depending on Eq and it is the energy 
that a particle of initial energy Eq has when it has reached depth N. 

An ion of energy E and charge z loses energy with depth at a r ate given by 
( Gouldl . 1972bl : Emslie, 1978 : Olive and Particle Data Grounl . 2014l l 


Here (3 is particle speed in units of the speed of light, K is given by 


(5) 


nieC^ 


Aeff{E,N) = xA+{l-x)A', (6) 

where x is the degree of ionisation and A and A' are the Coulomb logarithms 
appropriate to slowing on electro ns, and on neutral atoms respe c tively. For ions , 
expressions for A are given in e.g. Butler and Buckingham ( 1962 1: Gouldl (1972a h 
For A' we used the Bethe-Bloch form ( Olive and Particle Data Groun . 12014 '). 
neglecting the correction for the polarisation of the medium which is negligible 
in the low-density conditions of the solar atmosphere. We estimate the effects of 
the atmospheric chemical composition on energy loss by multiplying A' by 




a-H 


= 1.16 


where n sums over atomic species and A„, and a„ are the atomic number, 
atomic weight and abundance relative to hydrogen of species n. Also we gener¬ 
alise X = Ue/nH so that it is no longer quite the “degree of ionisation” and may 
take values > 1, dependent on the loca l electr o n dens ity. _ 

A and A' for electrons are given in ICould ( 1972al l ; lEmsliel ( 1978h . In using 
only Equation ([S]) we neglect pitch-angle scattering. Scattering in angle is unim¬ 
portant for fast ions and relativistic electrons - the only particles that can reach 
the LC7i laye r. For non-relativistic electrons it reduces the mean range by 2/3 
( Brown . 1972h . although dispersion means that electrons will stop and de posit 
energ y over a range of depths up to their full vertical stopping depth (IBail . 
1982l l. For the present purposes of estimation we neglect pitch-angle scattering 
completely. 

Figure [7] shows the value of Ap f /■ th r ougho ut the chromospheric portion of the 
C7 atmosphere of lAvrett and Loesei ( 2003 1. for protons of energies 1, 10 and 
100 MeV. The import ance of retaining the en e rgy-d ependence of A has been 
emphasised elsewhere ( MacKinnon and Tonerl . l2003ll . The variable degree of 
ionisation means that A^ff also depends strongly on position. Thus the energy 
loss rate of Equation ([S]) is not s e parab le in energy and depth, and analytical 
heating rates like those of lEmsli^ ( 1978ll may not be used (although an energy 
loss rate appropriate to a completely neutral atmosphere would give a good 
approximation below about 1000 km). 
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Figure 7. The effective Coulomb logarithm, through the atmosphere for protons of 

energy 1, 10 and 100 MeV. 


To evaluate the energy deposition 0 we need to assume a form for Fo{Eq), 
e.g. a single power law: 


^o(^'o) 




iS-1) 



(7) 


where Eq is the number of particles injected above energy if*. Then applying 
Equation JS]) in Equation (g]) gives 

0(iV) = 8.16xlO-^h^nH{N)EoEt-\S-l) H if-^(8 ) 

JErniniN) P [E) 

Again, E in the integrand of Equation ([5]) is understood to be a function of Eq 
and TV: the energy at depth of a particle that had energy Eq at injection. 
Since we cannot obtain analytical expressions for E{Eq,N) we calculate it via 
numerical integration of Equation with numerical interpolation as necessary 
between the tabulated points of a given atmosphere model, as needed for each 
of the abscissae of a numerical evaluation of the energy deposition (Equation 
dH)). More elaborate forms of Fo(ifo), e.g. the broken power-law deduced for 
electrons, are straightforwardly substituted in Equation (U) as necessary. 
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